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Abstract 

We employ data from precision electroweak tests and collider searches to derive constraints 
on the possibility that weak-singlet fermions mix with the ordinary Standard Model fermions. 
Our findings are presented within the context of a theory with weak-singlet partners for all 
ordinary fermions and theories in which only third-generation fermions mix with weak singlets. 
In addition, we indicate how our results can be applied more widely in theories containing exotic 
fermions. 
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1 Introduction 



The origins of electroweak and flavor symmetry breaking remain unknown. The Standard Model 
of particle physics describes both symmetry breakings in terms of the Higgs boson. Electroweak 
symmetry breaking occurs when the Higgs spontaneously acquires a non-zero vacuum expectation 
value; flavor symmetry breaking is implicit in the non-universal couplings of the Higgs to the 
fermions. However, the gauge hierarchy and triviality problems imply that the Standard Model 
is only an effective field theory, valid below some finite momentum cutoff. The true dynamics 
responsible for the origin of mass must therefore involve physics beyond the Standard Model. This 
raises the question of whether the two symmetry breakings might be driven by different mechanisms. 
Many theories of non-Standard physics invoke separate origins for electroweak and flavor symmetry 
breaking, and place flavor physics at higher energies in order to satisfy constraints from precision 
electroweak test and flavor-changing neutral currents. 

In this paper, we explore the possibility that flavor symmetry breaking and fermion masses 
may be connected with the presence of weak-singlet fermions mixing with the ordinary Standard 
Model fermions. Specifically, we consider theories in which some of the observed fermions' masses 
arise through a seesaw mechanism that results in the presence of two mass eigenstates for each 
affected flavor: a lighter mass eigensate whose left-handed component is predominantly weak- 
doublet, and a heavier one that is mostly weak-singlet. Such seesaw mass structures involving 
either third-generation fermions 0, [?J or all fermions || have played a prominent role in recent 
work on dynamical symmetry breaking. 

This work uses published experimental data to elicit constraints on the masses and mixing 
strengths of the exotic fermions. We both interpret our findings within the context of several 
specific models and indicate where our results can be applied more widely. Our initial approach 
is to study Z-pole and low-energy data for signs that the known fermions include a non-Standard, 
weak-singlet component. Previous limits of this type [|], ^, ^] have found that the mixing fraction 
sin 2 Omix can be at most a few percent for any given fermion species. As a complementary test 
we also look for evidence that new heavy fermions with a large weak-singlet component are being 
pair-produced in high-energy collider experiments. This can provide a direct lower bound on the 
mass of the new fermions. Most recent limits on production of new fermions focus on sequential 
fermions (LH doublets and RH singlets), mirror fermions (RH doublets and LH singlets), and vector 
fermions (LH and RH doublets) 0]. These limits need not apply directly to weak singlet fermions, 
as their production cross-sections and decay paths can differ significantly from those of the other 
types of fermions. 

We take as our benchmark a model || in which each ordinary fermion flavor mixes with a sep- 
arate weak-singlet fermion; this allows us to consider the diverse phenomenological consequences 
of the singlet partners for quarks and leptons of each generation. The low-energy spectrum is com- 
pletely specified, so that it is possible to calculate branching ratios and precision effects. Electroweak 
symmetry breaking is caused by a scalar, with flavor-symmetric couplings to the fermions. Fla- 
vor symmetry breaking arises from physics at higher scales that manifests itself at low energies in 
the form of soft symmetry-breaking mass terms linking ordinary and weak-singlet fermions. The 
fermions' chiral symmetries enforce a GIM mechanism and ensure that the flavor structure is pre- 
served under renormalization. Due to recent interest in using weak singlets to explain the mass 
of the top quark ||, we also analyze variants of our benchmark model in which only the third- 
generation fermions have weak-singlet partners. Furthermore, we indicate how our results can be 
applied to other theories with weak-singlet fermions. 

Since our benchmark model includes a scalar boson, it should be considered as the low-energy 
effective theory of a more complete dynamical model; specifically, at some finite energy scale, the 



I 



scalar <£, like the Higgs boson of the Standard Model, would reveal itself to be composite. That more 
complete model would be akin to dynamical top seesaw models f], 0, ||, which include composite 
scalars, formed by new strong interactions among quarks, and also have top and bottom quarks' 
masses created or enhanced by mixing with weak-singlet states. Those particular top seesaw models 
generally have multiple composite scalars when more than one fermion has a weak-singlet partner; 
these tend to be heavier than the single scalar in our models. Moreover, in the "flavor-universal" 
versions ||] generation-symmetry-breaking masses for the weak singlet fermions are the source of 
the differences between the masses of, say, the up and top quarks; the flavor structure of our models 
is different. Despite these differences, most of our phenomenological results are relevant to the top 
seesaw models. 

In the next section, we review the structure of our benchmark model, focusing on the masses, 
mixings, and couplings of the fermions. Section 3 discusses our fit to precision electroweak data |j| 
and the resulting general limits on the mixing angles between ordinary and weak-singlet fermions. 
We then use the constraints on mixing angles to find lower bounds on the masses of the new heavy 
fermion eigenstates. Section 5 discusses the new fermions' decay modes and extracts lower bounds 
on the fermion masses from LEP II [10, 11] and Tevatron |l^, [HJ data. Oblique corrections are 
discussed in section 6 and our conclusions are summarized in the final section. 



2 The Model 

At experimentally accessible energies, the models we consider have the gauge group of the Standard 
Model: SU(3)c x SU(2)w x U(l)y. The gauge eigenstate fermions include three generations of 
ordinary quarks and leptons, which are left-handed weak doublets and right-handed weak singlets 

tf> L = f^J , U R , D r U=(u,c,t), D = (d,s,b), 



D , 



L 



Ll = (Jj , Ir £={e,iM,r) . (2.1) 

In our general, benchmark model to each 'ordinary' charged fermion there corresponds a 'primed' 
weak-singlet fermion with the same electric charge^ 

U' L , R , D' l>r , i' LtR . (2.2) 

We will also discuss the phenomenology of more specialized models in which only third-generation 
fermions have 'primed' weak-singlet partners. 

The gauge symmetry allows bare mass terms for the weak-singlet fermions 

MuU' l U'r + M D D' L D' R + M e ? L l' R (2.3) 

and we take each of these mass matrices Mf to be proportional to the identity matrix. 
The model includes a scalar doublet field 



1 In principle, one could include weak singlet partners for the neutrinos as well. The neutrino phenomenology is 
largely separate from the issues treated here and will not be considered in this paper. 
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whose VEV breaks the electroweak symmetry. This scalar has Yukawa couplings that link left- 
handed ordinary to right-handed primed fermionic gauge eigenstates 



\u$l& U' r + \d^l$ D' r + \ e L L $ 4 . (2.5) 

The coupling matrices Xj are taken to be proportional to the identity matrix. The mass of the 
scalar is assumed to be small enough that the scalar's contributions will prevent unitarity violation 
in scattering of longitudinal weak vector bosons. 

Finally, there are mass terms connecting left-handed primed and right-handed ordinary fermions 

U' L muU R + D' L m D D R + e' L md R . (2.6) 

which break the fermions' flavor symmetries. We shall require the flavor-symmetry violation to be 
small: any mass my should be no greater than the corresponding mass Mf. This allows our model 
to incorporate the wide range of observed fermion masses without jeopardizing universality ||. 
As discussed in reference 101, this flavor structure is stable under renormalization. On the one 



hand, the flavor-symmetry-breaking mass terms (2.6) are dimension-three and cannot renormalize 



the flavor-symmetric dimension-four Yukawa terms ( |2.5| ). On the other, because all dimension- four 
terms (including the Yukawa couplings (|2.5|) ) respect the full set of global chiral symmetries, 

SU{^u, wDIr x SU(3)u, l x SU(3) d , l x SU(3) Ur x SU(3) Dr x 

SU(3) Ll j Ir x SU(3) Vl x SU(3) lR (2.7) 



they do not mix the mass terms ( |2.3| ) and (|2.q ) which break those symmetries differently. Fur- 
thermore, the global symmetries of this model lead to a viable pattern of inter-generational mixing 
among the fermions. Including the Mf terms ( |2.3| ) breaks the flavor symmetries to a form 

SU(S)i, L)WtD > x SU(3) Ur x SU(3) Dr x SU{3) Ll>v x SU(3) Ir (2.8) 

nearly identical to that of the Standard Model with massless fermions. Once the flavor-symmetry- 
breaking masses of equation (2.6) are added, the quarks' flavor symmetries are completely broken, 



leading to the presence of a CKM-type quark mixing matrix and an associated GIM mechanism that 
suppresses flavor-changing neutral currents. The lepton sector retains the U(l)'s corresponding to 
conservation of three separate lepton numbers. 



The ordinary and primed fermions mix to form mass eigenstates; for each type of charged 
fermion (/ = U, D, t) the mass matrix in the gauge basis is of the form 

This is diagonalized by performing separate rotations on the left-handed and right-handed fermion 
fields. The phenomenological issues we shall examine will depend almost exclusively on the mixing 
among the left-handed fermions. Hence, our discussion related to fermion mixing and its effects will 
focus on the left-handed fermion fields. For brevity, we omit "left" subscripts on the left-handed 
mixing angles and fields; we include "right" subscripts in the few instances where the right-handed 
mixings play a role. 

To evaluate the degree of mixing among the left-handed weak-doublet and weak-singlet fields, 
we diagonalize the mass-squared matrix (M'M). The rotation angle among left-handed fermions 
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is given by0 

B 2 

sm 2 cbf = l ; B = 2vXfM f (2.10) 

} 2A 2 + 2B 2 - 2AVA 2 + B 2 } 1 y ' 

A = M 2 + m 2 - v 2 X 2 

and the mass-squared eigenvalues are 

A± = I (m) + m) + v 2 X 2 f ) (l ± yl - (Av 2 X 2 m 2 )/{M 2 + m} + ^A 2 ) 2 ) . (2.11) 

Due to the matrix's seesaw structure, one mass eigenstate (f L ) has a relatively small mass, while 
the mass of the other eigenstate (f H ) is far larger. The lighter eigenstate, which has a left-handed 
component dominated by the ordinary weak-doublet state, 

/ i = cos0 / /-sin0 / /' , (2.12) 

corresponds to one of the fermions already observed by experiment. Its mass is approximately given 
by (for vXj < Mf and my < My) 

W ~^ + («A / )2 + m5 • (2 - 13) 
The heavier eigenstate, whose left-handed component is largely weak-singlet, 

/^ = sin < A / / + cos0 / / / (2.14) 

has a mass of order 

(mf) 2 M| + (vA/) 2 + m 2 , - (rn)) 2 . (2.15) 

The interactions of the mass eigenstates with the weak gauge bosons differ from those in the 
Standard Model because the primed fermions lack weak charge^]. The coupling of f L (f H ) to the 
W boson is proportional to cos 4>f (sin 0/) ; the right-handed states are purely weak-singlet and do 
not couple to the W boson. Thus the couplings of left-handed leptons to the W boson look like 
(since we neglect neutrino mixing) 



le 



smtiw 



(i L ~f^ cos </> e + I s sin W (2.16) 



When weak-singlet partners exist for all three generations of quarks, the left-handed quarks' cou- 
pling to the W bosons is of the form 

" (U^U^VudI °h W" (2.17) 



The 6x6 non- unitary matrix Vjj d is related to the underlying 3x3 unitary matrix Ajj p that mixes 
quarks of different generations 

y UD = ( Cu Aud Cd —Cjj A ud Sd \ ^ 2 
\ —Sjj Aud Cd Su Aud Sd ) 



2 To study the rotation angle among right-handed fermions, one diagonalizes (MM^) and obtains analogous results. 
3 For a general discussion of fermion mixing and gauge couplings in the presence of exotic fermions, see M. 
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through diagonal matrices of mixing factors 



Cjj = diag(cos <f) u , cos<p c , cos0 4 ) , Cjj = diag(cos (ftd, cos0 s , cos fa) , 

Sfj = diag(sm fa, sin0 c , sin^j) , Sd = diag(sin fa, sin</> s , sin 06 ) . 

The unitary mixing matrix A\jb, like the CKM matrix in the Standard Model, is characterized 
by three real angles and one CP-violating phase. But it is the elements of Vjjd which are directly 
accessible to experiment. While Vjjd is non-unitary, any two columns (or rows) are still orthogonal. 
The coupling of left-handed mass-eigenstate fermions to the Z boson is of the form 

ie (f L fH V/ ( cos 2 4> f T 3 -Q sin 2 6 W cos 0/ sin 0^ W f L \ 
smO w cos6 w KJ ,J )7fl \ cos^/sin^/Ts sin 2 0/T 3 - Q sin 2 9 W J \ f H J { > 

where T3 and Q are the weak and electromagnetic charges of the ordinary fermion. The right-handed 
states, being weak singlets, couple to the Z exactly as Standard Model right-handed fermions would. 
The scalar boson <I> couples to the mass-eigenstate fermions according to the Lagrangian term 

A/ (f L ,f H )ieft ( ~ C ° S t f Shl t f ' ri9M C ° S t f COS t f ™ ht )( t ) $ + h.c. (2.20) 
V - sm< ^/ sm <t>f,right sm0/ cos <j) fjHght j \ f H J V ' 

where (j) fright is the mixing angle for right-handed fermions. 

A few notes about neutral-current physics are in order. Flavor-conserving neutral-current decays 
of the heavy states into light ones are possible (e.g. [i H — > \i L v ^ . This affects the branching ratios 
in heavy fermion decays and will be important in discussing searches for the heavy states in Section 
5. Flavor-changing neutral (FCNC) processes are absent at tree-level and highly-suppressed at 
higher order in the benchmark model, due to the GIM mechanism mentioned earlier. For example, 
we have evaluated the fractional shift in the predicted value of T(b — > 57) by adapting the results 
in H|. As we shall see in sections 3 and 4, electroweak data already constrain the mixings between 
ordinary and singlet fermions to be small and the masses of the heavy up-type fermion eigenstates 
to be large (so that the Wilson coefficients 07(171/) that enter the calculation of T(b — > S7) are all 
in the high-mass asymptotic regime). The shift in T(6 — ► S7) is therefore at most a few percent, 



which is well within the 10% - 30% uncertainty of the Standard Model theoretical predictions [14] 
and experimental observations |l 



3 General limits on mixing angles 

Precision electroweak measurements constrain the degree to which the observed fermions can con- 
tain an admixture of weak-singlet exotic fermions. The mixing alters the couplings of the light 
fermions to the W and Z from their Standard Model values, as discussed above, and the shift in 
couplings alters the predicted values of many observables. Using the general approach of reference 
]l6fl , we have calculated how inclusion of mixing affects the electroweak observables listed in Table 
1. The resulting expressions for these leading (tree-level) alterations are given in the Appendix as 
functions of the mixing angles. We then performed a global fit to the electroweak precision data to 
constrain the mixing angles between singlet and ordinary fermions. The experimental values of the 
observables used in the fit and their predicted values in the Standard Model are listed in Table 1. 

To begin, we considered the benchmark scenario (called Case A, hereafter) in which all elec- 
trically charged fermions have weak-singlet partners ||. All of the electroweak observables given 
in Table 1 receive corrections from fermion mixings in this case. We performed a global fit for the 
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values of the 8 mixing angles of the fermions light enough to be produced at the Z-pole: the 3 
leptons, 3 down-type quarks and 2 up-type quarks. At 95% (90%) confidence level, we obtain the 
following upper bounds on the mixing angles: 

sin 2 (j) e < 0.0024 (0.0020) , sin 2 ^ < 0.0030 (0.0026) , sin 2 0r < 0.0030 (0.0025) 
sin 2 cf) d < 0.015 (0.013) , sin 2 <p s < 0.015 (0.011) , sin 2 cf) b < 0.0025 (0.0019) 
sin 2 (p u < 0.013 (0.011) , sin 2 <$> c < 0.020 (0.017) . (3.1) 

The 90% c.l. limits are included to allow comparison with the slightly weaker limits resulting from 
the similar analysis of earlier data in reference . 

The limits on the mixing angles are correlated to some degree. For example, most observables 
that are sensitive to d or s quark mixing depend on sin 2 (f)^ + sin 2 4> s . Indeed, repeating the global 
fit using the linear combinations (sin 2 4>d ± sin 2 </> s )/2 yields a slightly stronger limit for the sum 
(sin 2 (fid + sin 2 s )/2 < .0094 and a slightly looser one for the difference —0.0071 < (sin 2 — 
sin 2 s )/2 < .0195. This turns out not to affect our use of the mixing angles to set mass limits 
in the next section of the paper: limits on the d and s H masses arise from the more tightly- 
constrained 6-quark mixing factor sin 2 4>b instead. 

We, similarly, placed limits on the relevant mixing angles for three scenarios in which only third- 
generation fermions have weak-singlet partners. In Case B where the top quark, bottom quark and 
tau lcpton have partners, the 12 sensitive observables are TzjCr^, Rb, c ,e,ti,r, A b FB ,Ab, and R eT:)1T . 
The resulting 95% (90%) confidence level limits on the bottom and tau mixing angles are 

sin 2 <p T < 0.0018 (0.0014) , sin 2 <p b < 0.0013 (0.00084) (3.2) 

In Case C, where only the top and bottom quarks have partners, the nine affected quantities are 
Tz,Ch, Rb,c,e,tM,r, A b FB , and Ab- The sole constraint is 

sin 2 (j} b < 0.0013 (0.00084) (3.3) 

In Case D, where only the tau leptons have partners, only the six quantities Tz,(Jh, R T ,A FB , and 
R eT „r are sensitive, and the limit on the tau mixing angle is 

sin 2 (j) T < 0.0020 (0.0016) (3.4) 

These upper bounds on the mixing angles depend only on which fermions have weak partners, 
and not on other model-specific details. They apply broadly to theories in which the low-energy 
spectrum is that of the Standard Model plus weak-singlet fermions. 

4 From mixing angles to mass limits 

The constraints on the mixing between the ordinary and exotic fermions imply specific lower bounds 
on the masses of the heavy fermion mass eigenstates ( |2.15| ). We will extract mass limits from mixing 
angle limits first in the general case ||] in which all charged fermions have singlet partners, and 
then in scenarios where only the third generation fermions do. 

4.1 Case A: all generations mix with singlets 

Because the heavy fermion masses mj depend on v A/, Mf, and mj, we must determine the allowed 
values of all three of these quantities in order to find lower bounds on the m?. For the three fermions 
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Quantity 


Experiment 


SM 


Reference 


rv 

Zj 


2.4939 ± 0.0024 


2.4958 










41.491 ± 0.058 


41.473 








A,- (TV) 


0.1431 ± 0.0045 


0.1467 


IE 


1 




c \ T J 


0.1479 ± 0.0051 


0.1467 








At r 


0.1550 ± 0.0034 


0.1467 




1 (ref. 59) 




0.21656 ± 0.00074 


0.21590 








R c 


0.1735 ± 0.0044 


0.1722 








A b 


0.0990 ± 0.0021 


0.1028 








Ac 
FB 


0.0709 ± 0.0044 


0.0734 








Ah 


0.867 ± 0.035 


0.935 










0.647 ± 0.040 


0.668 








-72.41 ±.25 ± .80 


-73.12 ± .06 


0] (refs. 48, 49) 


Qw(Tl) 


-114.8 ± 1.2 ± 3.4 


-116.7 ± .1 


|j (refs. 48, 49) 




20.783 ± 0.052 


20.748 








20.789 ± 0.034 


20.748 






R T 


20.764 ± 0.045 


20.748 






A e 

FB 


0.0153 ± 0.0025 


0.01613 






AH 
FB 


0.0164 ± 0.0013 


0.01613 






AT 

FB 


0.0183 ± 0.0017 


0.01613 






As 
FB 


0.118 ± 0.018 


0.1031 ±.0009 






M W 


80.39 ± 0.06 


80.38 


HI 




q P v{v& — ► ve) 


-0.041± 0.015 


-0.0395±.0005 








5eA(^e -> i/e) 


-0.507 ± 0.014 


-0.5064±.0002 










0.3009±0.0028 


0.3040 ± .0003 










0.0328 ± 0.0030 


0.0300 


I 


1 




Rtt 


(1.230 ± .004)xl0~ 4 


(1.2352 ± .0005) xl0~ 4 




J ' 




Rr 


(1.347 ± .0082) xlO 6 


1.343 xlO 6 








RflT 


(1.312 ± .0087) xlO 6 


1.304 xlO 6 


I 


i\ 


1 



Table 1: Data used in fits to constrain mixing angles: experimentally measured electroweak ob- 
servables and their values within the Standard Model. 

of a given type, (e.g. e, /U, r), the values of Xf and Mf are common. The different values of nif 
arise from differences among the mj, and the form of equation ( p. 13 ) makes it clear that larger 
values of my correspond to larger values of mk. 

How can we ensure that the third-generation fermion in the set gets a large enough mass ? If 
we set mj to the largest possible value, ir\f = Mf, there is a minimum value of vXf required to 
make mf large enough. A smaller value of my would require a still larger value of vXf to arrive at 
the same wh. In other words, starting from ( |2,13j ), and recalling vXf < Mf we find 

vX f > \f2m L fz (4.1) 

where "/3" denotes the third-generation fermion of the same type as "/" (e.g. if "/" is the electron, 
"/3" is the tau lepton). The specific limits for the three types of charged fermions are: 

vX e > 2.5 GeV, vX D > 6.0 GeV vX v > 247 GeV (4.2) 
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Knowing this allows us to obtain a rough lower bound on the heavy fermion mass eigenstates. 
Since we require Mf >vXf and since the smallest possible value of mj is zero, we can immediately 
apply fly]) to ( fOED and find 

(mf f ^ A(m%) - (m)) 2 (4.3) 
For instance, the mass of the heavy top eigenstate must be at least 



m 



H > 



300GeV 



H - vom^«3UULieV (4.4) 

We can improve on these lower bounds in the following way. Because (mf) 2 is a monotonically 
increasing function of (vX) 2 , the minimum v\f, found above, yields the lowest possible value of 
mf . Thus, if we know what value of my should be used self-consistently with the smallest vXf, we 
can use (|2l5| ) to obtain a more stringent lower bound on rnf. The appropriate values 

(4.5) 



m 



/ = Mf (3rd generation) 



m, 



m) 



\ 



M 2 + v 2 X) 



v 2 X) 



(m L f ) 2 



(1st or 2nd generation) 



(4.6) 



follow from our previous discussion and from inverting equation ( |2.13| ), respectively. Because 
mf << Mf for the first and second generation fermions, our previous lower bound on mf for those 
generations is not appreciably altered. For the third generation we obtain the more restrictive 



(mf) 2 Z 



so that, for example, 



m 



H > 



V5; 



if ' ~ V u ^^t 

We can do still better by invoking our precision bounds on the mixing angles sin</>/. Recalling 
vXf < Mf and mj < Mf, allows us to approximate our expression (2.10) for the mixing angle as 



5(m/ 3 ) 2 



390 GeV 



(4.7) 



(4. 



sin (j)f 



vXfMf 



M 2 f + mj 



v 2 X 2 



(4.9) 



Further simplification of this relation depends on the generation to which fermion / belongs. For 
example, among the charged leptons, m e and are far smaller than Mi, while m r could conceivably 
be of the same order as Mg. Thus the limits on the leptons' mixing angles imply 



M e > Max 



vXp 



vXi vXp 



2 sin ( 



sin ( 



sin< 



(4.10) 



The strongest bound on Mi comes from sin(j) e ; that for Md, from sin that for Mjj, from sin< 



Mi > 



vX f 



smi 



> 51 GeV, 



M v > 



vXjj 



sm ( 



>2.2TeV, M D > 



vXp 
2 sin 6 



> 60 GeV 



(4.11) 



Combining those stricter lower limits on Mf with our bounds ( [4.1| ) on vXf and our expression 
for the heavy fermion mass ( (2.15|) gives us a lower bound on the mf for each fermion flavor. For 



1 

the third generation fermions we use ( |4.5| ) for the value of m f and obtain the 



> m^ 



m+ 



m* 



> m+ 



1 + 4/sin^ 
1 + 4/ sin 2 



b e > 73 GeV 
b u > 3.1 TeV 



> 



mg- Jl + l/sin 2 06 > 86 GeV 



c.l. lower bounds 

(4.12) 
(4.13) 
(4.14) 
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For the lighter fermions, we use equation (4.6) for the rr\f. Since mj << Mf in these cases, we find 



™?>™>jl ~ m^2 + 2/sin 2 e > 51GeV (4.15) 



mfJ2 + 2/sm 2 (t) u > 2.2 TeV (4.16) 



mg,m? ~ m^Y / 2 + l/2sin i ^ > 61 GeV . (4.17) 

The mass limits for the heavy leptons and down-type quarks are also represented graphically in 
figures |l] and ^. In figure [j], which deals with the leptons, the axes are the flavor- universal quantities 
Mi and v\j>. The shaded region indicates the experimentally allowed region of the parameter space. 



The lower edge of the allowed region is delimited by the lower bound on vXp of equation (4.2), as 
represented by the horizontal dotted line. The left-hand edge of the allowed region is demarked 
by the upper bound on the electron mixing factor, sin 2 <p e , as shown by the dashed curve with 
that label. The form of this curve, sin 2 4> e (Me, vXe) = 0.0024, was obtained numerically by using 
equation (2.11) for to put the unknown m e in terms of Mg, vX^ and the observed mass of the 



electron (m^ = .511 MeV) and inserting the result into equation ( |2.10| ). The curves for the muon 
and tau mixing angles were obtained similarly, but provide weaker limits on the parameter space 
(as shown by the dashed curves labeled sin 2 ^, and sin 2 </> T ). The lowest allowed values of the 
heavy fermion masses and are those whose curves intersect the tip of the allowed region; 
these are shown by the solid curves, obtained numerically by using equation (2.11) to replace the 



unknown m e , m T by the known m^,m^ in our expressions for and . Figure ^ shows the 
analogous limits on the mixing angles and heavy-eigenstate masses for the down-type quarks. 

We can also construct a plot of the allowed region of Mjj vs. vXjj parameter space. The lower 
edge comes from the lower bound on vXu and the left-hand edge, from the upper bound on sin^> u . 
We can then use the known value of m\ to calculate the size of the top quark mixing factor sin 2 fa 
at any given point in the allowed region. Numerical evaluation reveals 

sin 2 cf> t < 0.013 (0.011) (4.18) 

at 95% (90%) c.l. This is a limit on top quark mixing imposed by self-consistency of the model. 

In section |5j we will compare the mass limits just extracted from precision data with those 
derived from searches for direct production of new fermions at the LEP II and Tevatron colliders. 
The lower bounds on the masses of the heavy down-type quarks or charged leptons admit the 
possibility of those particles' being produced at current experiments. The heavy up-type quarks 
are too massive to be even singly produced at existing colliders. 

4.2 Cases B, C, and D: third-generation fermions mix with singlets 

If only third- generation fermions have weak-singlet partners, there are a few differences in the 
analysis that yields lower bound on heavy eigenstate masses. All follow from the fact that the 
lower bounds on the Mf (as in equation ( 4. 10| ) ) can no longer come from precision limits on the 



mixing angles of 1st or 2nd generation fermions (since those fermions no longer mix with weak 
singlets) . 

To obtain the precision bounds on the masses of b H and t h , we start by writing the lower limits 
on Mi and Md that come from the mixing angles: 

M,>-%, M D >£SL . (4.19) 
2 sin 0t 2 sin (by. 
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Figure 1: A graphical representation of the coupling and mass limits for the e H and t h states in the Mi 
vs. vXi plane. The shaded region indicates the experimentally allowed region of the parameter space. The 
quantities on the axes are flavor-blind; separate curves for different lepton flavors are shown. Each dashed 
curve shows the upper bound on sin 2 (f>g for one lepton and excludes the region to its left. Each solid curve 
shows the lower bound on mf for one lepton (the limits for e and ijl are identical, as discussed in the text 
and equation(4.1£ ). The horizontal dotted line indicates the minimum allowed value of v\? (4.1). 
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Figure 2: A graphical representation of the mass limits for the d H and b H states in the Mo vs. vAp plane. 
The shaded region indicates the experimentally allowed region of the parameter space The quantities along 
the axes are flavor-blind; separate curves for different quark flavors are shown. Each dashed curve shows 
the upper bound on sin 2 (f>o for one quark and excludes the region to its left. Each solid curve shows the 
lower bound on for one quark (the limits for d and s arc identical, as discussed in the text and equation 



( 4.17 ). The horizontal dotted line indicates the minimum allowed value of v\d (4.1) 
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The factor of 2 in the denominator arises because the mixing angles belong to a third-generation 
fermion (so that mj = Mf). We therefore find 



mf > + l/sin 2 ^ T (4.20) 

mf > m^A/l + 1/sin 2 (4.21) 



In Case B where all third-generation fermions mix with weak-singlet fermions, the mixing angle 
limits ([T^) based on the twelve sensitive observables yield 95% c.l. lower bounds 



mf > 42 GeV (4.22) 
mf > 119 GeV (4.23) 



In Case C, where only third-generation quarks have partners, (|3.3| ) which was obtained by a fit to 
the nine aaffected observables, gives 

mf > 119 GeV (4.24) 

while in Case D, where only tau leptons have partners, (|3.4j ) based on six affected precision elec- 
troweak quantities implies 

mf > 40 GeV (4.25) 

Compared with the limits in Case A, we see that the lower bound on mf is strengthened because 
the precision limit (|3.2|, |3.3j ) on sin 2 ^ is more stringent. In contrast, the lower bound on mf 
is weakened because the bound now depends on a third-generation instead of a first-generation 
mixing angle: equation ( 4.20) ) is approximately mf > m^/sm.(p T whereas equation ( |4.12| ) was 



roughly mf > 2m^ / sin <j) e . 

Note that in theories where the top is the only up-type quark to have a weak-singlet partner, 
such as Cases B and C, the only bound on mf comes from equation ( |4.8| ). While this is far weaker 
than the limit in Case A, it still ensures that the heavy top eigenstate is too massive to have been 
seen in existing collider experiments, even if singly produced. 



5 Limits on direct production of singlet fermions 



While interpreting the general mixing angle limits in terms of mass limits requires specifying an 
underlying model structure, it is also possible to set more general mass limits by considering 
searches for direct production of the new fermions. The LEP experiments have published limits 
on new sequential charged leptons |lc|| |ll|; the Tevatron experiments have done the same for new 



quarks [12|[13|. In this section, we adapt the limits to apply to scenarios in which the new fermions 
are weak singlets rather than sequential. 



5.1 Decay rates of heavy fermions 

A heavy fermion decays preferentially to a light fermionF] plus a Z, W, or $ boson which subsequently 
decays to a fermion-antifermion pair (see figure ||) [] 

4 Even where a heavy fermion is kinematically allowed to decay to another heavy fermion, the rate is doubly- 
suppressed by small mixing factors (sin^/) and, consequentially, negligible. 

5 In this section we confine our analysis to relatively lig ht scalars, with mass below 130 GeV. For heavier scalar 
one should include the scalar decays to W and Z pairs |20| and the resulting 5-fermion final states of heavy fermion 
decays. We expect this to yield only a small change in the results of our quark-sector analysis and essentially no 
alteration in our results for heavy lepton decays, due to the large kinematic suppression when m$ >> mf ~ Mw ■ 
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Figure 3: Scalar ($) and weak boson (V = Z or W) decay modes of a heavy fermion (f H ). 



At tree-level, and neglecting final state light fermion masses, we obtain the following partial 
rates for vector boson decay modes of the heavy fermions 
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fj fk fl 
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(5.1) 
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where V represents a Z or W boson, while Ty and My are, respectively, the vector boson's decay 
rate and mass. Function F(x,y) is presented in appendix ||. The vertex factors cfj (c^) are, as 



shown in figure ||, the fl* fj'V (f^ff'V) couplings which may be read from equations (2.16) - ( |2.19| ). 

Our results for the charged-current decay mode agree with those presented in integral form in 
1 21]. Moreover, equation yields the standard asymptotic behaviors in the limit of heavy fermion 
masses far above or far below the electroweak bosons' masses (see appendix |B|). Since some of our 



heavy fermions can, instead, have masses of order 80-90 GeV, we use the full result (5.1) in our 
evaluation of branching fractions and search potentials. 

Similarly, we find the partial rate for the scalar decay mode to be 
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where T$ and M$ are the decay rate and the mass of the scalar boson, $. Function G(x, y) and 
additional details are given in appendix [B]. The vertex factors c* (c^) are, as indicated in figure 
||, the fl 1 fj$ (/^/; L $) couplings which may be read off of equation fl2.20p . 

We have numerically evaluated the couplings of the light fermions to the scalar^, Z, and W as 
functions of the Mf and the v\f. In the region of the model parameter space that is allowed by 
precision electroweak measurements, we find that these couplings are within 1% of their Standard 
Model values. Therefore, in this section of the paper, we approximate the $// and Vff couplings 
for the light fermions by the Standard Model values. This allows us to express our results for 

6 To evaluate the m ixing among right-handed ferm ions which appears in the fermion-scalar couplings, we derive a 
relation analogous to (2.11) and apply equation 2.13 so that sin 2 4>j, right is written in terms of known light fermion 
masses, the Mf and the v\f. 
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branching fractions and searches in the simple Mf vs. vXf planes for the up, down, and charged- 
lepton sectors. In this approximation, the recent LEP lower bound on the mass of the Higgs boson 
p2|| , Mh > 95.3 GeV, applies directly to the mass of the $ scalar in our model: 

M # > 95.3 GeV (5.3) 

The branching ratios for the decays of the heavy leptons are effectively flavor-universal, i.e. the 
same for e H , fJL H , and t h . The charged-current decay mode dominates; decays by Z emission are 
roughly half as frequent and decays by $ emission contribute negligibly for mf < M$ . In the limit 
where the heavy lepton masses mf are much larger than any boson mass, the branching ratios 
for decays to W, Z, and approach 60.5%, 30.5%, and 9%, respectively. The branching fractions 
for heavy lepton decays are shown in figure § as a function of heavy lepton mass mf , with M$ 
fixed at 100 GeV and v\i set equal to IvrJ^. As the branching ratios have little dependence on the 
small mixing factors sin^y (as we argue in more detail in the following subsection), they are also 
insensitive to the value of vXf. 

The branching fractions for decays of the heavy down-type quarks display a significant flavor- 
dependence. Those for the d H and s H are essentially identical and resemble the branching fractions 
for the heavy leptons. However, charged-current decays of b H with a mass less than 255 GeV (the 
threshold for decay to an on-shell top and W) are doubly Cabbibo-suppressed, so that the b H 
branching ratios do not resemble those of the other down-type quarks. Generally speaking, a b H 
of relatively low mass decays almost exclusively by the process b H — ► Zb L . For m b H larger than 
255 GeV, the decay b H — > t L W dominates and the Z-mode branching fraction is only about half 
as large. If mf is above M$ + m^ but below 255 GeV, the scalar decay mode becomes significant 
(in agreement with reference p3]| ). If the scalar mass lies above 255 GeV, the scalar decay mode is 
much less important. In the asymptotic regime, where is much greater than mt or any boson's 
mass, the branching ratios for decays to W, Z, and $ approach 49%, 25%, and 26%, respectively. 
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5.2 Heavy leptons at LEP II 

The LEP II experiments have searched for evidence of new sequential leptons, working under the 
assumptions that the new neutral lepton N is heavier than its charged partner L and that L decays 
only via charged-current mixing with a Standard Model lepton (i.e. B(L — > U(W*) = 100%). 



Recent limits from the OPAL experiment at yfs = 172 GeV [10] and from the DELPHI experiment 



at y/s = 183 GeV [11] each set a lower bound of order 80 GeV on the mass of a sequential charged 
lepton. 

To illustrate how the LEP limits may be applied to our weak-singlet fermions, we review OPAL's 
analysis. The OPAL experiment searched for pair-produced charged sequential leptons undergoing 
charged-current decay: 

e+e" -> L + L- -» u^W+W (5.4) 

Their cuts selected final states in which at least one of the W* bosons decayed hadronically. Events 
with no isolated lepton were required to have at least 4 jets and substantial missing transverse 
momentum; those with one or more isolated leptons were required to have at least 3 jets, less 
than 100 GeV of visible energy, and substantial missing transverse momentum. The efficiencies for 
selecting signal events were estimated at 20-25% by Monte Carlo. With 1 candidate event in the 
data set and the expectation of 3 Standard Model background events, OPAL excluded, at 95% c.l., 
sequential leptons of mass less than 80.2 GeV, as these would have contributed least 3 signal events 
to the data. 

The heavy leptons in the models we are studying have different weak quantum numbers than 
those OPAL sought. This alters both the production rate and the decay paths of the leptons. 
The production rate of the £ H should be larger than that for the sequential leptons. The pure 
QED contribution is the same, as the heavy leptons have standard electric charges; the weak- 
electromagnetic interference term is enhanced since the coupling to the Z is roughly sin 2 9\y > 



rather than sin &w — 0.5 < as in the Standard Model. By adapting the results of reference [24] 



to include the couplings appropriate to our model, we have evaluated the production cross-section 
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Figure 6: Production cross-section for a heavy lepton that is mostly weak-singlet as a function of lepton 
mass and mixing angles. Each family of curves represents one value of y/s (solid = 172 GeV, dashed = 189 
GeV, dotted = 192 GeV). The separate curves within each family show the effect of changing the value of 
the small mixing angles (top curve: sin 2 <j)e = 0.1 for all leptons, middle curve: sin 2 fa = for all leptons, 
bottom curve: sin 2 fa, = 0, sin 2 4>^ T = 0.1). 



for heavy leptons at LEP II. Our results are shown in figure |6| as a function of heavy lepton mass 
for several values of \fs and lepton mixing angle. 

On the other hand, the likelihood that our heavy leptons decay to final states visible to OPAL 
is less than it would be for heavy sequential leptons. In events where both of the produced l H 
decay via charged-currents, about 90% of the subsequent (standard) decays of the W bosons lead 
to the final states OPAL sought - just as would be true for sequential leptons. But the heavy 
leptons in our model are not limited to charged-current decays. In events where one or both of 
the produced £ H decay through neutral currents, the result need not be a final state visible to 
OPAL. If there is one W and one Z in the intermediate state, about 36% of the events should 
yield final states with sufficient jets, isolated leptons and missing energy to pass the OPAL cuts. 
At the other extreme, if both £ H decay by <3? emission, there will be virtually no final states with 
sufficient missing energy, since <3? decays mostly to bb. The other decay patterns lie in between; for 
intermediate ZZ ($Z, §W) we expect 28% (19%, 30%) of the events to be visible to OPAL. The 
total fraction of pair-produced heavy leptons that yield appropriate final states is the sum of these 
various possibilities: 

Bdecay = 0.896 {B w f + 0.280 {B z f + 2 ( 0.361 B w ■ B z 

+ 0.190 B z + 0.306 B$ -BwO (5-5) 

where Byy, Bz, and B$ are the heavy lepton branching fractions for the W, Z and scalar decay 
modes respectively, as calculated in section 5.1 (and shown in Figure 4). 

In models (cases B and D) where there is only one species of heavy lepton (r ), setting a mass 
limit is straightforward. We note that 

(^production ' Bdecay — ^events / £ ' C> (^ - ^) 
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Figure 7: Lower bounds on heavy lepton mass in models with one flavor of heavy lepton (Cases B and D). 
The heavy dashed curve encloses the region of the plot allowed by the combination of electroweak and direct 
search data. The limits on sin</> T 2 (dashed lines) and v\i = (dot-dash line) are as in Figure 1. Solid 

lines are contours of constant heavy lepton mass m^; dotted curves are contours of constant numbers of 
OPAL signal events. The lower bound the lepton mass comes from the overlap of the limiting N — 3 dotted 
line with the — 79.8 GeV solid line. In calculating branching fractions, the scalar mass was set to 100 
GeV. 



where, as in OPAL's analysis, the integrated luminosity is C = 10.3 pb^ 1 , the signal detection 
efficiency^] is e 20%, and the number of (unseen) signal events is N events ~ 3. Thus an upper 
bound on the number of signal events implies an upper bound on a pr oduction • B decay Inserting 
the branching fraction for £ H £ H pairs to visible final states, B 'decay: as in equation ( |5.5| ) yields an 
upper bound on the production cross-section. Since we have already calculated the cross-section 
{<fproductiori{\fs = 172 GeV)) as a function of heavy lepton mass, we can convert the bound on 
^production into a a 95% c.l. lower bound on 



H. 



777 



f > 79.8 GeV. (5.7) 



This is a great improvement over the bounds of order 40 GeV, ( |4.22| ) and (4.25), we obtained 



earlier from precision electroweak data in cases B and D where the tau is the only lepton to have 
a weak-singlet partner. 

Our new lower bound on further constrains the allowed region of the Mi vs. vXe parameter 
space, as illustrated in figure [?]. Contours on which the heavy tau mass takes on the values = 
70, 79.8 and 92 GeV are shown as a reference and to indicate how a tighter mass bound would 
affect the size of the allowed region. 

In case A, where e, n, and r all have singlet partners, the contributions from all three heavy 
leptons to the signal have to be taken into account. While the e H and \i H have nearly identical 

7 Our use of OPAL's 20% signal efficiency is conservative. OPAL considered pair-production of sequential leptons 
that decay via charged currents. About one-tenth of the time, both W's decay leptonically; these ttvvvv final states 
would be rejected by OPAL's cuts. In considering cases where one or both of our heavy leptons decay via neutral 
currents, we have not included the analogous llvvvv events. Thus a higher fraction of the events we did include 
should pass OPAL's cuts. 
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Figure 8: Lower bounds on heavy lepton mass in models with three flavors of heavy lepton (Case A). The 
heavy dashed curve encloses the region of the plot allowed by the combination of electroweak and direct search 
data. The limits on sin0 T 2 (dashed lines) and vXg = (dot-dash line) are as in Figure 1. Solid lines are 

contours of constant heavy lepton mass - those for are vertical, those for are curved. Dotted lines 
are contours of constant number of OPAL signal events; the cusp shows where the pair-production threshold 
for t h is crossed. The mass limit for e H , /i H comes from the overlap of the N events = 3 dotted line with 
the rn^ = 84.9 GeV solid line; at these values of vXe and Mg, the t h is above OPAL's pair-production 
threshold. In calculating branching fractions, the scalar mass was set to 100 GeV. 



masses and decays, the t h has slightly different properties. By adding the contributions from all 
three flavors of heavy lepton, drawing the contour corresponding to N even t s = 3 on the Mivs.vXi 
parameter space, and comparing this with contours of constant mf for each species, we obtain the 
95% c.l. lower bounds on all three heavy lepton masses, as shown in figure |8] 

mf tli > 84.9 GeV (5.8) 

mf > 93.9 GeV . (5.9) 

Note that the bound on comes simply from internal consistency of the model (the values of v\g 
and Mi are flavor- universal) , since it lies above OPAL's pair-production threshold. These bounds 
are a significant improvement over those we obtained from precision data, i.e. ( [4.15 ) and Q4.12 ). 



While calculating the lower limits on the mf required us to assume a value for M$ (to evaluate 
Bdecay), the result is insensitive to the precise value chosen. As noted in section 5.1, in the allowed 
region of the v\i vs. Mi plane, LEP's lower bound on the Higgs boson's mass applies to $ so that 
mm.{mf) < min(M$). In this case, B{i H — ► $>£ L ) is negligible. 

Our limits are also insensitive to the precise values of the small lepton mixing angles sin<^. The 
production rate has little dependence on sinr/)^ because the £ H £ H Z coupling ( [2.19|) is dominated by 
the "— Qsin 2 ^" term. What little dependence there is on smcp£ decreases as 2m¥ approaches \/s, 
and the mass limits tend to be set quite close to the production threshold. Moreover, the branching 
fractions for the vector boson decays of the £ H have only a weak dependence on sin^. Both the 
charged- and neutral-current decay rates are proportional to sin 2 <p£ (and the rate for decay via 
Higgs emission is negligible), so that the mixing angle dependence in the branching ratio comes 
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only through factors of cos 2 4>i which are nearly equal to 1. As a result, our lower bounds on the 
heavy fermion masses will stand even if improved electroweak measurements tighten constraints on 
the mixing angles. 

Because the mass limit tracks the pair-production threshold, stronger mass limits can be set by 
data taken at higher center-of-mass energies. Figure ^ shows CF proa <uction as a function of the heavy 
lepton mass for several values of -^/i and sin 2 <f>£. As data from higher energies provides a new, more 
stringent upper bound on (^production ' E 'decay, one can read an improved lower bound on the heavy 
lepton mass from figure [f| 

More generally, one can infer a lower mass limit on a heavy mostly-weak-singlet lepton from 
other models using the same data by inserting the appropriate factor of Bd ecay in equation ( |5.6[ ), 
For models in which the mixing angles between ordinary and singlet leptons are small and in which 
B(£ H — > <&£ L ) is small, our results apply directly. This would be true, for example, of some of the 
heavy leptons in the flavor-universal top seesaw models ||. 

Since the lower bound the LEP II data sets on the mass of the heavy leptons is close to the 
kinematic threshold for pair production, it seems prudent to investigate whether single production 

e + e- -» £ H £ L (5.10) 

would give a stronger bound. Single production proceeds only through Z exchange (the jf H f L 
coupling is zero). Moreover, equation ( 2,19| ) shows that the Z£ H £ L coupling is suppressed by a 



factor of sin<^; given the existing upper bounds on the mixing angles ( |3. I| )-( p74| ) , the suppression 
is by a factor of at least 10. As a result, only a fraction of a single-production event is predicted to 
have occurred (let alone have been detected) in the 10 pb _1 of data each LEP detector has collected 
- too little for setting a limit. 

5.3 Heavy quarks at the Tevatron 

New quarks decaying via mixing to an ordinary quark plus a heavy boson would contribute to 
the dilepton events used by the Tevatron experiments to measure the top quark production cross- 



section |12|][13]. We will use the results of the existing top quark analysis and see what additional 
physics is excluded. If evidence of new heavy fermions emerges in a future experiment, it will be 
necessary to do a combined analysis that includes both the top quark and the new fermions and 
that examines multiple decay channels. 

Here, we use the dilepton events observed at Run I to set limits on direct production of new 
largely- weak-singlet quarks (our q H ). These new quarks are color triplets and would be produced 
with the same cross-section as sequential quarks of identical mass. However, their weak-singlet com- 
ponent would allow the new states to decay via neutral-currents as well as charged-currents. This 
affects the branching fraction of the produced quarks into the final states to which the experimental 
search is sensitive. 

The D0 and CDF experiments searched for top quark events in the reaction 

pp^QQ^ qWqW -» qqlv^'ve (5.11) 

by selecting the final states with dileptons, missing energy, and at least two jets. Di-electron and di- 
muon events in which the dilepton invariant mass was close to the Z mass were rejected in order to 
reduce Drell-Yan background. The top quark was assumed to have essentially 100% branching ratio 
to an ordinary quark (q) plus a W, as in the Standard Model. The D0 (CDF) experiment observed 
5 (9) dilepton events, as compared with 1.4 ± 0.4 (2.4 ± 0.5) events expected from Standard model 
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backgrounds and 4.1 ± 0.7 (4.4 ± 0.6) events expected from top quark production. Thus, D0 
(CDF) measured the top production cross- section to be 5.5 ± 1.8 pb (8.2±|;| pb). 

In using this data to provide limits on the production of heavy quarks in our models, we consider 
dilepton events arising from top quark decays to be part of the background. Hence, from D0 (CDF), 
we have 5 (9) dilepton events as compared with a background of 5.5 ± 0.8 (6.8 db 0.8) events. At 
95% confidence level, this implies an upper limit of 5.8 (9.6) on the number of additional events 
that could have been present due to production and decays of new heavy quarks. 

How many q H would be produced and seen ? The q H have the same QCD production cross- 
section as a Standard Model quark of the same mass. The q H can decay by the same route as the 
top quark ( |5.11| ). About 10% of the charged-current decays of pair-produced q H would yield final 
states to which the FNAL dilepton searches were sensitive. The neutral current decays of the q H 
reduce the charged-current branching fraction B(q H — > q L W), but will not, themselves, contribute 
significantly^] to the dilepton sample since dileptons from Z decays are specifically rejected and the 
couplings to e and fj, are extremely small. Then we estimate the fraction of heavy quark pair 
events that would contribute to the dilepton sample as 

Bdecay = (Bw) 2 (5-12) 

where Bu is the fraction of W pairs in which both bosons decay leptonically and Bw = B(q H — » 
q L W) is calculated in section 5.1 and shown in figure 5. 

The number of dilepton events expected in a heavy-quark production experiment with luminos- 
ity C and detection efficiency for dilepton events e is 

N qH = a qH ■ C ■ e ■ B decay (5.13) 

Similarly in top searches the total number of events is 

N l = a l -C-e-Bcc- B u (5.14) 

where B cc is the fraction of top quark pairs decaying via charged currents. 

In comparing the number of events expected for produced top quarks with those for q H pairs, the 



values of e and C are the same; furthermore, B cc of equation ( 5.14 ) is essentially 100% . Therefore 
we may write 

a q " (B w f = a t ^ r (5.15) 

Using the values which the CDF and D0 experiments have determined for the three quantities on 
the right-hand side (cf. previous discussion), we find 

a qH (B w ) 2 < 7.8 pb (DO) 

< 12.0 pb (CDF). (5.16) 

The dilepton sample at the Tevatron is sensitive only to the presence of the d H or s H quarks in our 
models. The u H , c H and t H are, according to equations (|4.8|), (|4.13| ) and ( 4.16 ), too heavy to be 



produced, while the b H decay dominantly by neutral instead of charged currents, due to Cabibbo 
suppression. Hence, this search tests only the models of case A, in which the light ordinary fermions 
have weak-singlet partners. 



A Higg s-like scalar with a mass of order 130-150 GeV could have a relatively large branching fraction to two W 
bosons pel . This might allow some neutral-current decays of q 11 to contribute to the dilepton sample and change 
our mass bounds slightly. 
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Figure 9: A graphical representation of the Tevatron limits on heavy D quark masses. The region allowed 
by precision electroweak tests (cf. fig 2) is bounded from above by the dashed line (sin^f, 2 ) and from below 
by the dot-dashed line (vXd)- Dotted lines represent curves of constant <J prod ■ (B w ) 2 as in equation (5.15). 
Solid vertical lines are curves of constant rn ds (m^); their overlap with the dotted lines defines the lower 
bound on the heavy fermion masses. We set M$ = 100 GeV in calculating branching fractions. 



Since the pair-production cross-section for q H is the same as that for a heavy ordinary quark, 
we use the cross-section plots of reference [25] and our calculated branching fraction B^ ecay (5.12) 



to translate equation ( 5.16| ) into lower bounds on heavy fermion masses. For Case A, in which both 



the d H and s H quarks can contribute to the dilepton sample, we find (with M<j> = 100 GeV): 



> 


153 GeV 
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161 GeV 


(CDF) 



(5.17) 

which are significantly stronger than those obtained from low-energy data in section 4 and also 



stronger than the published limits on a fourth- generation sequential quark [26]. Note that since the 
b H decays almost exclusively via neutral-currents due to Cabbibo suppression of the charged-current 
mode, the lower bound on is, once again, an indirect limit implied by internal consistency of 
the model. In the scenarios where only third-generation fermions have weak partners (Cases B and 
C), we can obtain no limit on mff . 

More generally, one can use the same data to infer an upper limit on the pair-production cross- 
section for heavy mostly-weak-singlet quarks from other models by inserting the appropriate factor 
of B\v in equation (|5.16[) . After taking into account the number of heavy quarks contributing, one 
can use the cross-section vs. mass plots of p5j] to determine lower bounds for the heavy quark 



masses. For example, our cross-section limits ( 5.16 ) apply directly to the heavy mostly-singlet 



quarks in the dynamical top-seesaw models that are kinematically unable to decay to scalars and 
decay primarily by charged-currents. The corresponding mass limit depends on how many such 
quarks are in the model. 
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Figure 10: Comparing data on oblique corrections to theoretical predictions. Relative to the reference 
[m t = 173.9 GeV, m H = 300 GeV, a~ 1 (M z ) — 128.9], the cross shows the best experimental fit to S and 
T; the solid ellipses are at the corresponding 68%, 90% and 95% confidence levels |29j for two degrees of 
freedom. The labels on both the solid and dotted ellipses indicate A\ 2 relative to the experimental best-fit 
point (cross). The heavy dotted curve shows how the predicted value of S and T in the Standard Model 
varies as the scalar mass to$ is varied by steps of 10 GeV (see text); lower masses are to the left. The value 
of m$ corresponding to the lowest x 2 (smallest dotted ellipse) is 80 GeV. 



6 Oblique Corrections 

The presence of new singlet fermions present in our models will shift the S and T parameters 
from their Standard Model values. In this section, we evaluate these changes and explore the limits 
they impose on the fermion masses and couplings and the mass of the scalar, <£. This analysis of 
one-loop oblique corrections turns out to complement the analysis of tree-level effects on precision 
data performed in section 3: the oblique corrections most strongly limit the top quark mixing angle 
which the earlier analysis could not directly constrain. 

In calculating the values of S and T predicted by our models, we started from the results of 



[29], which cite the experimental values of S and T relative to the reference point [rat = 173.9 
GeV, rriH = 300 GeV, a~ x {Mz) = 128.9]. We included the appropriately weighted variations of 
mt and a" 1 and obtained the minimal combined x 2 field on the S — S re f vs T — T re j plane; we 
simultaneously obtained the corresponding mt(S,T) and a -1 (S,T) that minimize x 2 f° r each pair 
of S and T parameters. The minimal combined x 2 is presented in in figure [l(]; the solid ellipses 
represent joint 68.3%, 90%, and 95.4% c.l. limits on S and T with variations in mt and a -1 included. 



Next, within the Standard Model we allowed the Higgs mass to vary [28] from 40 GeV to 1 TeV in 
steps of 10 GeV and obtained the "best fit Higgs curve" shown in figure [l(]; the circled points are 
at 100, 200, 300, 400, and 500 GeV (smaller masses to the left). The dotted ellipses in the figure 
are contours of constant minimal combined x 2 whose intersections with the "best fit Higgs curve" 
define the best fit value and 68.3%, 90%, and 95.4% c.l. limits on Higgs mass. These values are 
respectively 80 GeV (in good agreement with f§), 190 GeV, 310 GeV, and 400 GeV. 

We then added the effects of the extra fermions on S and T. The contribution of the singlet 
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fermions to S was calculated numerically using the formalism described in 
to T was found analytically 



The contribution 



, 32] by summing the vacuum-polarization diagrams containing the 
heavy and light mass-eigenstate fermions present in the model of interest. For example, in models 
containing weak-singlet partners for only the t and b quarks, we find that the contribution of the 
t H , t L , b H and b L states to the T parameter is (in agreement with J?]]) 
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(6.1) 



where Tg is the Higgs contribution, and cj (sj) is an abbreviation for cos (ft f (sin eft j). To isolate 
the extra contribution caused by the presence of the weak-singlet partners for the t and b quarks, 
we must subtract off the amount which t and b contribute in the Standard Model [32]: 
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(6.2) 



Note that ( |6.1| ) correctly reduces to ( |6.2j ) in the limit where singlet and ordinary fermions do not 
mix (sin 2 (ftt , sin 2 (ft b — > 0). From the form of equation fl6.1| ), we see that experimental bounds on 
the magnitude of T will constrain relatively heavy extra fermions to have small mixing angles. 

To illustrate how oblique effects constrain non-standard fermions, we begin by including a weak- 
singlet partner only for the top quark; that is, we send sin 2 (ft b — » in equation (6.1). For a given 



scalar mass m$, we add to the Standard Model S and T, the additional contribution caused by 
mixing of an ordinary and weak-singlet top quark. For the T parameter, this extra contribution is 
the difference between expressions ( |6.lD and (|6.2| ) with sin 2 cft b = 0. By construction, for s 2 — > the 
new contributions to the S and T parameters both go to zero (i.e. 5S = 5T = 0). When mixing is 
present (sf ^ 0), one has 5S < and ST > 0, and the predicted values of S and T lie above the 
"best fit Higgs curve" 

We deem "allowed" the values of m^ and sin 2 (ftt for which the final values of S and T fall inside 
the dotted ellipse labeled A% 2 = 5.25 - the 90% c.l. ellipse for the Standard Model alone. In other 
words, we require that the model including new physics agree with experiment at least as well as 
the Standard Model. This allows us to trace out a region of allowed heavy top mass and mixing 
for different values of m$ , as illustrated in figure |ll|. Note that the presence of non-zero mixing of 
ordinary and singlet top quarks enables a heavier scalar to be consistent with the data[]. 

As a complementary limit on m^ and sin 2 <ftt , we note that the discussion in section 4 requires 



771 



H 



> mfJl + 1/sin^ 



(6.3) 



For a discussion of related issues for the Standard Model Higgs boson see pgj ]. 
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Figure 11: Lower bound on heavy top mass mf as a function of heavy top mixing sin 2 <j) t . Based on the 
oblique corrections, for m$ = 100 GeV, the mass and mixing must fall below the solid curve; for m$ = 350 
GeV, they must fall in the band between the dashed curves; for m$ = 520 GeV, they must lie within the 
dotted curve. The additional lower bound on mf 6.3 is represented by the heavy solid curve; the allowed 
region is to the right of this curve, leading to the constraints (|6.4[) and 



That is, for a given amount of mixing, the heavy top mass must lie above some minimum value. 



Combining these limits yields the allowed region in the mixing vs. mass space in figure 11. For 
example, 



For rriA = lOOGeV, 



sin 



For m = 350GeV, 



sin 



> 



< 
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1450GeV 
.015 

1040GeV 
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(6.4) 
(6.5) 



As illustrated in figure [ll], if the scalar 's mass, M$, rises above 520 GeV, the regions of top mass 
and mixing allowed by oblique corrections by equation ( |6.3[ ) cease to intersect; this provides an 
upper bound on the scalar mass. 

To apply oblique-correction constraints to our models, we need to include weak-singlet partners 
for quarks other than the top quark. Since these fermions contribute little to S p7|], we can illustrate 
the effects of including other singlet fermions by showing how they affect the T parameter. First, 
we include the singlet partner for the b quark, as in equation (|6.1|). We can interpret the result 



using figure 12, which shows the value of T within the coupling- mass plane for the up-sector quarks. 
For reference, dotted nearly-vertical curves of constant heavy top mass are shown. The main 
contents of the figure are the three sets of curves labeled ST = [0.3, 0.1, 0], where 5T is the 
contribution due to mixing between ordinary and singlet fermions. Within each set, the separate 
curves correspond to different values of the heavy b mass and mixing. The solid curve obtains for 
5 TeV and sin 2 <f> h = 0.00090; the dashed curve, for = 5 TeV and sin 2 fa = 0.00040; the 



m 



doted curve, for = 0.55 TeV and sin 2 (f>b = 0.00027. Looking at the region where is of order 
a few TeV, we see that the influence of the b-quark is small. Including the effects of partners for 
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Figure 12: The effect of b-quark mixing on T. Contours of constant <5T are shown in the up-sector 
Muy.s.vXu parameter space. Three representative values (ST = 0, 0.1, 0.3) are shown for several values of 



heavy b-quark mass and mixing. Solid curves correspond to = 5TeV and sin 2 



0.00090 (equivalently, 



M D = 5TeV, v\ D = 150GeV). Dashed curves are mf = 5TeV and sin 2 4> b = 0.00040 {M D = 5TeV, 
v\ D = lOOGcV). Dotted curves are for mf = 0.55TeV and sin 2 <j> h = 0.00027 {M D = 0.5TeV, v\ D = lOGeV). 
The nearly-vertical dotted curves of constant m^ 1 are shown for reference. 



the other fermions yields a generalized version of equation 6.1 and similar results. Thus the lower 
bounds on we found earlier by considering only mixing for the top quark will not be much 
altered by including mixing for the other quarks, as in our models A, B, and C. 



7 Conclusions 



Precision electroweak data constrains the mixing between the ordinary standard model fermions 
and new weak-singlet states to be small; our global fit to current data provides upper bounds on 
those mixing angles. Even when the mixing angles are small, it is possible for most of the exotic 
mass eigenstates which are largely weak-singlets to be light enough to be accessible to collider 
searches for new fermions. We have analyzed in detail a class of models in which flavor-symmetry 
breaking is conveyed to the ordinary fermions by soft symmetry-breaking mass terms connecting 
them to new weak-singlet fermions; such models have a natural GIM mechanism and a flavor 
structure that is stable under renormalization. By calculating the branching rates for the decays 
of the heavy mass-eigenstates (which are significantly influenced by their being primarily weak- 
singlet in nature) we have been able to adapt results from searches for new sequential fermions to 
further constrain our models. We find that direct searches at LEP II now imply that the heavy 
leptons £ H must have masses in excess of 80-90 GeV; those limits are not sensitive to the precise 
values of the small mixing angles. Current Tevatron data indicates that heavy quark states d H 
and s H could be as light as about 140-150 GeV, while the mostly-weak-singlet b H must weigh at 
least 160-170 GeV. In addition, the new fermions' contributions to the oblique corrections allow 
the scalar $ to have a relatively large mass (up to about 500 GeV) while remaining consistent 
with the data. Oblique corrections also constrain the mixing and mass of the the heavy top state 
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which is mostly weak-singlet; in particular, must be at least 1 TeV. Finally, we have indicated 
how our phenomenological results may be generalized to related models, including the dynamical 
top-seesaw theories. 



A Appendix: Mixing effects on electroweak observables 

This appendix contains the expressions for the leading-order (in mixing angles) changes to elec- 
troweak observables in the presence of fermion mixing. The expressions were derived using equations 



(2.16 - 2.19) and the general approach of reference 16] 



Ar z /r|- 



SM 



0.603 (sin 2 fa + sin 2 fa) - 0.072 sin 2 fa 

0.3535 (sin 2 fa + sin 2 (j) s + sin 2 fa) - 0.287 (sin 2 fa + sin 2 fa) 



(A.l) 



Acr h/V h 



-1.409 sin 2 fa + 0.736 sin 2 c/> M + 0.072 sin 2 fa 

0.1515 (sin 2 fa + sin 2 <j) s + sin 2 (f> b ) - 0.124 (sin 2 fa + sin 2 . 



(A.2) 



AA T (P T ) = AA e (P T ) = AA LR = -0.5180 sin 2 e + 1.2870 sin 2 ^ (A.3) 



AR h /R% 



SM 



-0.0295 (sin 2 <f> e + sin 2 ^) + 0.505 (sin 2 <p d + sin 2 <j> a ) 
1.78 sin 2 (j) b + 0.411 (sin 2 </) u + sin 2 (j> c ) 



(A.4) 



AR C /R, 



SM 



= 0.0605 (sin 2 4> e + sin 2 + 0.505 (sin 2 <f> d + sin 2 4> s + sin 2 <f> h ) 
+ 0.411 sin 2 4> u - 1.999 sin 2 4> c 



(A.5) 



AA b FB = -0.3300 sin 2 4> e + 0.8500 sin 2 ^ - 0.0161 sin 2 4> h 



AA C FB = -0.1785 sin 2 e + 0.6665 sin 2 fa - 0.0875 sin 2 , 



AA b = 0.1052 (sin 2 <p e + sin 2 fa) - 0.1472 sin 2 fa 



AA C = 0.5719 (sin 2 fa + sin 2 fa) - 0.7997 sin 2 



(A.6) 
(A.7) 
(A.8) 
(A.9) 



AQ w (Cs) = 72.7663 sin 2 fa - 0.7239 sin 2 fa + 211.0024 sin 2 fa - 187.9988 sin 2 fa (A.10) 



AQ w {Tl) = 111.396 sin 2 fa - 4.920 sin 2 fa + 327 sin 2 fa - 285 sin 2 fa (A.ll) 
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AR e /Rf M = 2.275 sin 2 e + 0.130 sin 2 0^- 0.505 (sin 2 (f) d + sin 2 <p s + sin 2 <f> b ) 

- 0.411 (sin 2 (j> u + sin 2 C ) (A.12) 

AR^/R^ M = 0.130 sin 2 <p e + 2.275 sin 2 <^ — 0.505 (sin 2 <^ + sin 2 <p s + sin 2 ^) 

— 0.411 (sin 2 + sin 2 <f> c ) (A.13) 

AR T /Rf M = 0.130 (sin 2 </> e + sin 2 ^) + 2.145 sin 2 (f) T (A.14) 
- 0.505 (sin 2 <j) d + sin 2 S + sin 2 4> b ) - 0.411 (sin 2 (j) u + sin 2 C ) 

AA e FB = -0.1230 sin 2 (f> e + 0.3070 sin 2 (p^ (A. 15) 

AAp B = 0.0920 (sin 2 </> e + sin 2 </> M ) (A. 16) 

= 0.0920 sin 2 (f> e + 0.3070 sin 2 </> M - 0.2150 sin 2 (f> T (A. 17) 

AA^ B = -0.3300 sin 2 cf> e + 0.8500 sin 2 ^ - 0.0161 sin 2 <f> 3 (A.18) 

AMhz/M^ = 0.1065 (sin 2 e + sin 2 M ) (A. 19) 

A 5e y(zve -► i/e) = 0.1720 sin 2 </> e - 0.3650 sin 2 M (A.20) 

A5 e A(^e -► ve) = 0.5000 sin 2 </> e - 0.5060 sin 2 ^ (A.21) 

Ag\(vN -► i/A) = 0.1220 sin 2 e + 0.7260 sin 2 ^ - 0.4280 sin 2 d - 0.3445 sin 2 4> u (A.22) 

Ag 2 R (vN -► i/A) = -0.0425 sin 2 e + 0.0179 sin 2 M (A.23) 

AR n /R% M = - sin 2 e + sin 2 where R n = ffi — ^ (A.24) 

r(vr -> /xz^J 

AR eT /R? M = sin 2 0„ - sin 2 T where # er = ^ ez ^) ( A2 5) 

AR^/R™ = sin 2 e - sin 2 T where i? MT = ^ — (A.26) 
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B Appendix: Details of heavy fermion decays 



This appendix contains details relevant to the heavy fermion decays discussed in section 5.1. 

At tree-level and neglecting final-state light fermion masses, the kinematic factors F(x,y) and 
G(x,y) referred to in the text have the following form: 



F(x,y) = — {2 x(2-x) + [3(x-l)+ y 2 ]A(x,y) 
x A 

+ [(x-l) 2 (x + 2) + 3y 2 (x-2)]B(x,y)} 



(B.l) 



G(x,y) 



: {x (4x - 3) + [x (4 
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(B.2) 



where A and B are given by 
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(B.3) 
(B.4) 



To check our general expressions for the decay rates, we evaluated their behavior in the limiting 
cases where the decaying heavy fermion is either much more massive or much less massive than the 



vector or scalar boson involved in its decay. Equations 5.1 and B.l for vector-boson decays yield 
asymptotic behavior 
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where V may be either Z or W. Equations 5.2 and B.2 for scalar-boson decays yield 
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